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1. Introduction 

A tiling of a checkerboard with dominoes is a way of putting dominoes on the board so 
that no square of the board is uncovered and no two dominoes overlap. Given a local 
pattern (see figure [l] for examples), a location in the board, and the shape and size of 
the board, how many tilings of the board have the given pattern at the given location? 
(Alternatively, we can substitute "bond" for "domino" and "particle" for "square", and 
ask for the probability of local patterns in a system of particles each of which bonds with 
exactly one of its neighbors.) 

Suppose that the squares of the board are very small compared to the board itself. For 
some board shapes, the probability of finding a pattern at a given location will be the same 
for almost all locations. This is the case for the square board. (See figures || to @, where 
tiles are colored according to their direction and parity for the sake of clarity; see figure 
for the coloring scheme.) There are some boards, however, for which the probability 
does depend on the location. Consider, for example, the Aztec diamond, that is, the board 
whose boundary is a square tilted 45 degrees (figures and ||) . In random tilings of the 
Aztec diamond, we usually find brick-wall patterns outside the inscribed circles, and more 
complicated behavior inside the circle. (See figures || to O.) 

The probabilities of local patterns in a rectangular board were computed recently |(| . Until 
now, there was no other board for which the probabilities of all local patterns were known. 
Many experiments and some important partial results (lj had shown that, as already stated, 
the probabilities of patterns in the Aztec diamond depend on location. This qualitative 
difference between the Aztec diamond and the rectangular board made the former as worthy 
of analysis as the latter. The main result of this work is an expression for the probability of 
any local pattern in a random tiling of the Aztec diamond. The expression is a determinant 
of size proportional to the number of squares in the pattern, just like Kenyon's expression 
|6| for the probabilities in the rectangural board, 

Main Result 1. The probability of a pattern covering white squares v\, 1)%, ■ • ■ Vk and black 
squares wi, 11)2, ■ ■ • Wk of an Aztec diamond of order n is equal to the absolute value of 

H^^j)\ l]=l2 ... k . 

The coupling function c{v, w) at white square v and black square w is 

Xi-l 

2 " ^2 kr (- 7 ' n ' Vi ~ ^ kr (y'i -l,n-l,n- (j + xl i - aft)) 

for xli > Xi and 

n 

-2~ n ^2 kr C?' 71 > Vi - !) kr (y'i - 1, - 1, n - (j + xU - Xi)) 

for xli < Xi, where (xi,yi) and (a;/,-,^) are the coordinates of v and w, respectively, in the 
coordinate system in figure |7J, and the Krawtchouk polynomial kr(et, b, c) is the coefficient 
ofx a in {l-x) c ■ {l + x) b ~ c . 

Our line of attack is as follows. 

1. Reduce the problem of finding probabilities of patterns to an enumerative problem; 




m 



(a) 



(b) 



Figure 1 . A few examples of local patterns 



(C) 




Figure 2. These two patterns have the same probability of being found 
in a random pattern at any given place 




FIGURE 3. Random tiling of 
square of side 20 




Figure 4. Random tiling of 
square of side 40 



2. Reduce the enumerative problem to a simpler one involving Aztec diamonds with two 
holes rather than arbitrary even-area holes; 

3. Compute the weighted number of tilings of an Aztec diamond with two holes. 



The first two steps involve known techniques, and were already considered to be a plausible 
strategy by other researchers. The third step is new. 



Henry Cohn is currently analyzing the case of the board with infinitely small squares by 
approximating the sum of Krawtchouk polynomials in our main result as an integral for 
n — > oo. His results will be presented in a later, joint version of this paper. 



Figure 5. Random tiling of 
square of side 80 



horizontal, 

white square on the left 

(when board colored as a checkerboard) 



horizontal 

black square on the right 



vertical, 

white square on the top 



■ vertical, 
black square on the top 

Figure 6. Shading chart 



2. The Kasteleyn Matrix 

What is the probability of finding a pattern in a tiling of a given board? It is equal, by 
definition, to the number of tilings of the board with the given pattern, divided by the total 
number of tilings of the board. Clearly, the number of tilings of the board with the given 
pattern depends only on the squares occupied by the pattern, and not by how they are tiled 
by the pattern: there is a one-to-one correspondence between tilings with the given pattern, 
and tilings of the board with the squares covered by the pattern removed. (See figure 0.) 
Thus, what we want to know is the number of tilings of the board with the squares covered 
by a pattern removed, divided by the number of tilings of the board. 



Figure 7. Aztec diamond of 
order 4, as a board 




Figure 8. Aztec diamond of 
order 4, as a graph 




Figure 9. Random tiling 

of Aztec diamond of order Figure 10. Random tiling 

(side) 20 of Aztec diamond of order 40 



Figure 11. Random tiling 
of Aztec diamond of order 60 



Figure 12. Random tiling 
of Aztec diamond of order 80 




Figure 13. Random tiling of Aztec diamond of order 160 

Kasteleyn || showed that the number of tilings of a board can be expressed as the absolute 
value of a determinant with half as many rows as there are squares in the board. Thus, 
in a sense, our problem is already solved: since determinants can be computed in time 
polynomial on the number of rows, we can compute the probability of any pattern in a board 
in time polynomial on the size of the board. Unfortunately, there are two disadvantages to 
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this approach. The first one is that, if the Aztec diamond in question has side n, we will 

2 

have to compute a determinant of side \ , and that demands a considerable amount of time 
and space. Even more seriously, a sequence of determinants of varying size is very hard to 
analyze asymptotically. We would like to analyze such a sequence in order to determine 
the probability of a pattern near a point for Aztec diamonds of very high order, that is, of 
very fine "grain" (see figure |l3| ). Thus, a Kasteleyn determinant is not good enough. A 
determinant whose size depended only on the size of the pattern, and not on the size of the 
board, would be much easier to manipulate. 

In this section, we will prove that the number of tilings of a board equals the absolute 
value of the determinant of its Kasteleyn matrix, and then show how this implies that the 
probability of a pattern in a random tiling of a board is equal to a minor of the inverse 
of the Kasteleyn matrix of the board. This minor has side proportional to the number of 
squares in the pattern. Finally, we will show that the problem of finding an entry in the 
inverse of a Kasteleyn matrix can be reduced to an enumerative problem concerning an 
Aztec diamond with one black hole and one white hole. 

We will henceforth refer, not to boards, but to their dual graphs, which can be seen as 
having a vertex at the center of every square and an edge perpendicular to every edge 
between two squares. For example, the Aztec diamond will mean for us the object in figure 
H and not the object in figure This convention will simplify graph-theoretical arguments 
considerably. 

The results in this section are in part a modern formulation of Kasteleyn's work, and in 
part a codification of local folklore, as crafted and passed down by R. Kenyon, J. Propp, D. 
Wilson and others. 

2.1. The Kasteleyn- Wilson matrix. Consider a finite subgraph G of the infinite square 
lattice with as many white as black vertices, where the infinite square lattice is colored as 
a checkerboard. Let (i>\, ■ ■ ■ ,v n ) be its white vertices and (wi, ■ ■ ■ ,w n ) its black vertices. 
(Any ordering from 1 to n can be chosen.) The Kasteleyn- Wilson matrix 

K((vi, ■ ■ ■ ,v n ), (wi, ■ ■ ■ ,w n )) 

(or K(G), by abuse of language) is defined to be |aij|™, where a,-j is 

• 0, if {vi,Wj} is not an edge of G; 

• 1, if {vi,Wj} is a horizontal edge of G; 

• (— 1) , if {vi,Wj} is a vertical edge going from row I to row I + 1 and there are k 
vertices in row I to the left of the edge. 

We will show that the number of perfect matchings of G is equal to the absolute value of the 
determinant of K{G). As could be exprected, the Kasteleyn- Wilson matrix is only one of 
several Kasteleyn matrices K{G) 1 that is, determinants whose absolute values are equal to 
the perfect matchings of G. We will prove the enumerative property of the Kasteleyn- Wilson 
matrix because it is true for any subgraph G of the infinite square lattice, and not only for 
the Aztec diamond. Fortunately, the proof can be applied to other Kasteleyn matrices with 
minimal cases. In fact, in the last subsection of this section, and in following sections, we 
will use the following convention, which is more convenient for our purposes than Wilson's: 
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• 0, if {v i} Wj} is not an edge of G; 

• 1, if {vi,Wj} is a horizontal edge of G; 

• (— l) fc , if {v i} Wj} is a vertical edge going from row I to row I + 1 and there are k 
vertical edges from row I to row I + 1 to the left of the edge. 

This convention is valid for the Aztec diamond and for any other subgraph G of the infinite 
square lattice such that any lattice vertex inside any loop in G is also a vertex of G. 
Thus, only the material in the following two subsections is valid for any subgraph of the 
infinite square lattice. The rest of this work is specific to the Aztec diamond, although the 
techniques used are applicable to other boards. 

2.2. Why does K give us the number of perfect matchings? We want to show that 
the number of perfect matchings of a subgraph G of the infinite square lattice is equal to 
the absolute value of the determinant of K(G). 

We can express det(K) as 

n 

sgn(7r) JJa ii7r(i) , (1) 

7r€P({l,2,...,n}) *=1 

where P ({1, 2, n}) is the set of all permutations of {1, 2, • • • n). Define a map / from the 
set of all perfect matchings of G to P({1, 2, n}) as follows. Any perfect matching can be 
expressed in the form 

{{vi,w k(1) },{v 2 ,w k(2) },- ■ ■ ,{v n ,w k{n) }}, (2) 

where k is a map from {1, 2, • • -n} to itself. The map / takes {{v\, w k (i)}, • • • , {v n , w fe („)}} 
to k. It is clear that A; is a permutation; otherwise there would be unpaired vertices, as well 
as vertices belonging to more than one pair. Thus, / is well defined. Moreover, 

1. / is injective: If two matchings had the same map k, they would be the same matching. 

2. every element k of the image of / satisfies (JIILi a i,k(i)) 7^ 0: if Vi,w k (i) is an edge, 
then a,i.k(i) must be non-zero. 

3. if a permutation 7r of 1,2, ...,n satisfies Jl"=i a i,ir(i) 7^ 0, then, for every 1 < i < n, 
(vi,w n (i)) is a valid edge. Moreover, for i\ ^ i 2 , ^ 7i"(*2), and thus (v^ , w n ^) 
and (vi 2 , W7r(i 2 )) do not have any vertices in common. Therefore 

{{vi,w k{1) },{v 2 ,w k{2) }, - ■ ■ ,{v n ,w k(n) }} (3) 

is a perfect matching. 

Hence / is a one-to-one and onto map from the set of all perfect matchings of G to the set 
of all permutations k of 1, 2, n satisfying fllLi a i.k(i) 0- Therefore there are as many 
perfect matchings as there are non-zero terms in 

n 

det(K)= ^2 sgn(7r) JJa i)ff (i). (4) 

7reP({l,2,-,n}) i=l 

Every non-zero term is equal to either 1 or —1. To prove that the absolute value of det(_R') 
equals the number of perfect matchings, we have to show that all non-zero terms have the 
same sign. 
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Let M = {{v i: uvj)}}-^! and Ml = {{v i7 be two perfect matchings of G. By the 

definition of perfect matching, every vertex of G is in one edge of M and in one edge of Ml. 
It follows that every vertex of G is either in two edges or in no edges of (MUM/) — (Mn Ml). 
Therefore (MUM/) — (MflM/) consists entirely of loops, that is, it is a collection of disjoint 
sets of the form 

{{«i a , u>ji }, {wj! ,v i2 }, {v l2 ,w h },--- {v lm , w jm }, {w im , % }}. (5) 

(See figure O.) If a vertex is in two edges of (M U M/) — (M fl M/), one of these two edges 
must be in M, and the other one in M/. We can assume without loss of generality that 
{vi 1 ,Wj 1 } is in M, and hence {wj t ,Vi 2 } is in Ml, {vi 2 ,Wj 2 } is in M, and so on. Then, on 
one hand, ji = 7r(£j) for 1 < I < m, and, on the other hand, ji = 7r/(i; + i) for 1 < I < m — 1, 

im = 7r/(«i). Hence i 2 = ((^O^ 1 ° 7r )(*i)> *3 = (C^) -1 ° t)(«2), • ■ ■ ii = ((t^) -1 ° ^(v))- 
Thus every loop in (MUM/) — (MflM/) induces a cycle in (7r/) _1 o7r. It is easy to see that, 
conversely, for every cycle in (7r/) _1 o tt there is a loop in (M U Ml) — (M n Ml) that induces 
it. Because the sign of a permutation is equal to the product over all its cycles of (—1) to 
the power of the length of the cycle minus one, and because a cycle has length equal to half 
the number of edges of the loop inducing it, we have 

sgn((^)- 1 o 7 r) = n(-l) lonW/2 " 1 , (6) 



where L is the set of all loops of (M U Ml) — (M fl Ml) and len(^) is the number of edges 
in loop £. From this equation, from 

and from the fact that all 1 or (—1), it follows that the result we want to prove in 

this section, namely, 

n n 

sgn(7r) Y[ a»,,r(i) = sgn(7i7) a ii7r , (i) , (8) 



is equivalent to 



Now, 



JJ ( _ 1) ic„ W/2 -i = Y[a lMl) • n^T/W- (9) 

i£L i=l i=l 



H ffl»,ir(t) " JJ a i,7r/(i) = JJ ( H a i,7r(i) ' H a»,*7(i)) (10) 

»=i i=i feL ie/(^) 

, where I(£) is the set of indices i of all white vertices in loop £. Therefore it is enough for 
us to prove that 

(-l) le "W/2-l = . -Q fl . ^ (n) 

iei(e) iei(e) 
for every loop £ in (M U Ml) - (M n M/). 

,7r(i) ' I\i£i(e) a i,-!T'(i) is the product a^j over all 1 < i,j < n such that {vi,Wj} 
is an edge of loop £. Since ajj = 1 for {wi,^} horizontal, we can restrict the product 
to 1 < i,j < n such that {vi,Wj} is a vertical edge of loop I. What is, specifically, the 
product of <2ij over all 1 < i,j < n such that {vi,Wj} is a vertical edge of loop £ having 
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Figure 14. The union of two perfect matchings of the Aztec diamond, 
minus their intersection 



a vertex on row y and another in row y a + 1, where yo is a given? Let {(an, j/o), ( x i,Uo + 
1)}, {(x2,yo), (x 2 ,yo + 1)}, • • • {(x m ,y ), (x m ,y + 1)} be all such edges. (We are referring 
vertices by their Cartesian coordinates.) Since the horizontal line {(t, y + 1) : t e (— oo, oo)} 
crosses the loop an even number of times, m must be even. Let m — 2 ■ m . Since all a^j 
are 1 or (—1), we have 

m m a 

TT u Jo(x2j,yo), Ji(x 2 j, 2/0+1) / 10 x 

a Jo(a; : ,-,yo),Jl(^,yo + l) - 11 " i U^J 

j- =1 j = i "Jo(a:23-i,J/o),Ji(a!2 3 '-i,!/o+l) 

where Jo(z,w) is the index of the white vertex of coordinates (z,w), and J\(z,w) is the 
index of the black vertex of coordinates (z, w). Since aj ( Xi ,yo),Ji(xi,yo+i) ls equal to the 
number of vertices of G on row yo and to the left of Xi, and aj (x i+ i,j/ ),^i(x» + i,i/o+i) * s ec l ua l 
to the number of vertices of G on row yo and to the left of Xj+i, 

a ./o(3:i + i,yo),Ji(a;i + i,yo + l) 
a J (x i ,j/o),Ji(a: i ,yo + l) 

is equal to the number of vertices of G on row y and to the left of x i+ i but not of Xi. This 
is the same as the number of vertices of the infinite square grid on row y and to the left 
of x i+ i but not of Xi, minus the number of vertices in the grid but not in G, on row y 
and to the left of x i+ i but not of a^. The number of vertices of the grid to the left of x i+ i 
but not of Xi is equal to the number of squares of the grid contained between the edges 
((xi,y),(xi,y + l)) and {{x i+1 , y), (x i+1 , y + 1)). 

It is clear that, when i is even, the squares between the edges ((xi, y), (xi, y +1)) and 
((xi + i,y), (xi + i,y + 1)) arc in the interior of the loop, as are the vertices on row y of the 
grid which do not belong to G and which are the left of x i+1 but not of Xi. (Since any 
vertex on the loop, and, specifically, Xi and a;,+i, must be in G, we can hencefort refer to 
these vertices as "the vertices on row y of the grid which do not belong to G and which 
lie between Xi and x i+ i" .) Conversely, every square between rows y and y + 1 and in the 
interior of the loop lies between edges {{xi,y), (xi,y + 1)) and ((xi+i,y), (xi+i,y + 1)) for 
some odd i, and, moreover, every vertex on row y of the grid, not in G, and in the interior 
of the loop lies between Xi and Xi+i for some odd i. Hence the number of squares between 
edges ((xi, y), (xi, y + 1)) and ((a;,+i, y), (xi+i, y + 1)) for i odd equals the number of squares 
which are both between rows y and y + 1 and in the interior of the loop, and, furthermore, 
the number of vertices not in G lying on row y between Xi and Xi+i for i odd equals the 
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number of vertices not in G lying on row y and in the interior of the loop. Thus 

m mo 

a Jo(^,i/ ),Ji(^,yo+l) - 11 " " I 14 ) 

j—l j =1 U Jo(x2j-l,yo),Jl(x2j-l:yo + l) 

equals (—1) to the power of the number of squares which are both in the interior of the loop 
and between rows y and y + 1 minus the number of such vertices in the interior on row y. 
Hence, the product of this expression over all rows, that is, 



Hie/ a i,Tr/(i) 



(15) 



equals (—1) to the power of the number of squares in the interior of the loop minus the 
number of vertices lying on row y and in the interior, but not in G. 

Therefore 

^■-fa-.n^Mo (16) 

is equal to (—1) to the power of the length of the loop, divided by 2, minus 1, plus the 
number of squares inside the loop, minus the number of vertices in the grid and inside the 
loop, but not in G. Pick's theorem states that, for polygons whose vertices belong to a 
square grid, A = I + B/2 — 1, where A is the area enclosed by the polygon, / is the number 
of grid points inside the polygon, and B is the number of grid points on the boundary of 
the polygon. Hence our expression is equal to (—1) to the power of the number of vertices 
inside the loop, minus the number of vertices in the grid and inside the loop, but not in G. 
This is the same as the number of vertices of G inside the loop. Now, the vertices inside the 
loop are matched only among themselves, not with the vertices outside the loop, in either 
M or Ml. For a set of vertices to be matched only among themselves, there must be an 
even number of vertices in the set. Therefore the number of vertices of G inside the loop 
must be even, and (—1) raised to the power of this number must be 1. 

Hence 

( _i)^-in^ffl°^> =1 (17) 

1 !*€/(<) a h*<(i) 

as we desired to prove. We conclude, by taking the product over all loops I £ L, that 



sgn(7i7) n i= i Oi,7r/(i) 

Therefore the terms of 

n 

det(K)= sgn{n)l[a iMi) (19) 

7r£P({l,2,...,n}) i=l 

all have the same sign. It follows that the number of perfect matchings, that is, of permu- 
tations 7r for which n"=i a i.^(i) 1S non-zero, is equal to the absolute value of det(K). 
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2.3. Why does K^ 1 give us the probabilities of patterns? Let us have a graph G 
and a subgraph H whose every vertex has degree 1. There is a one-to-one correspondence 
between perfect matchings of G having H as a subgraph and perfect matchings of G — H. 
(By G — H we mean the graph (V, E), where V is the set of all vertices of G not in H and 
E is the set of all edges of G between two vertices of G not in H.) Therefore the probability 
that a random perfect matching of G have the set of edges of if as a subset is equal to the 
number of perfect matchings of G — H divided by the number of perfect matchings of G. 

What is, then, the number of perfect matchings of G — HI One way to compute it is 
to construct a Kasteleyn- Wilson matrix for G — H. Another way, which will soon prove 
its virtues, is to take the minor of the Kasteleyn- Wilson matrix for G resulting from the 
deletion of the rows and columns corresponding to H. Certainly, this minor is not the same 
as the Kasteleyn- Wilson matrix for G — H; the signs of the matrix entries are different. 
Nevertheless, the absolute value of the determinant of the minor is equal to the number 
of perfect matchings of G — H. To prove this, we need to do the same as in the previous 
section, namely show that all non-zero terms of the expression of the determinant as a 
sum over permutations have the same sign. If we proceed in the same way as before, we 
arrive at a point where the only difference is that we have (—1) to the power of the number 
of vertices of G inside a loop instead of the number of vertices of G — H inside the same 
loop. This difference is no difference if there is an even number of vertices of H inside the 
loop. Since the vertices of H inside the loop cannot be connected with the vertices of H 
outside the loop, the vertices of H inside the loop are paired among themselves. Thus we 
have, as we wanted, that there is an even number of vertices of H inside the loop. This is 
enough for us to show that the terms in the expression of the minor of K{G) as a sum over 
permutations do not cancel. 

Therefore the number of perfect matchings of G — H is equal to the absolute value of 
the determinant of the minor of K(G) lacking the rows and columns corresponding to 
the vertices of H. If we assign the same probability to every perfect matching of G, the 
probability that a random perfect matching will have H as a subgraph will be equal to 
the number of perfect matchings of G — H divided by the number of perfect matchings of 
G, that is, to the absolute value of the determinant of the aforementioned minor of K{G) 
divided by the absolute value of the determinant of K(G). By Jacobi's rule, a corollary 
of Cramer's rule, this is equal to the absolute value of the determinant of the minor of 
(K (G) _1 ) T consisting of those rows and columns omitted from the minor in the numerator, 
that is, of rows 1 < b\ < b 2 < ■ ■ ■ b m < n and columns 1 < c\ < c 2 < • • • c m < n, where 
Vbi , Vb 2 , ■ ■ ■ vt m and w Cl , w C2 , • • • w Cm are the vertices of H. 

This is a clear improvement over the expression ^pj^jjjp- Instead of dealing with deter- 
minants of the size of G, we deal with a determinant of the size of H. As we explained in 
the introduction, we are interested in finding the probability of small subgraphs, or "local 
patterns" , in a large graph G. We want to find what happens when we have an infinite se- 
quence of G's whose number of vertices goes to infinity. Now it is enough for us to examine 
a fixed number of entries in each (K{Gi)~ 1 ) T , and determine their asymptotic behavior as 
i — > oo. 

2.4. How can we find the entries of K~ x for an Aztec diamond by counting 
perfect matchings? Suppose that we have to compute the determinant of the minor of 
K(G)~ 1 consisting of rows 1 < b\ < b 2 < ■ ■ ■ < b m < n and columns 1 < c\ < c 2 < • • • < 
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c m < n. To do so, we have to compute the entry {{K{G)~ 1 ) T )b i , Cj for 1 < i, j < m. In other 
words, in order to compute the probability of any local pattern, it suffices to be able to 
compute an arbitrary entry of the inverse Kasteleyn matrix. If we have a sequence of graphs 
{Cfe}fc°=i (such as, for example, Aztec diamonds of higher and higher order) we will know 
the asymptotic behavior of the probabilities of local patterns if we know the asymptotic 
behavior of the entries in the sequence of matrices {K^Gk) -1 }^^- 

We can, of course, compute a first minoij^of K(G), and then apply Cramer's rule, whenever 
we want an entry of (K(G)~ 1 ) T . Unfortunately, it seems very hard to obtain asymptotic 
expressions directly from the minors. We will reduce the problem of computing the minor 
of K(G) resulting from the deletion of row i and column j to an enumerative problem whose 
solution we will be able to represent in a form other than a determinant. It would seem, 
at first sight, that we can use the same kind of argument we used to answer the previous 
question, and show that such a minor is equal (up to sign) to the number of all perfect 
matchings of G with white vertex i and black vertex j deleted. Unfortunately, this is not 
the case. A first minor without row i and column j is equal to a sum whose number of 
terms is equal to number of perfect matchings of G with white vertex i and black vertex 
j deleted. The problem is that, while every term has absolute value 1, not every term 
has the same sign. It is easy to show, by the same kind of reasoning we employed in our 
answer to the first question, that the matter of whether or not two terms have the same 
sign can be determined by examining the loops in the superimposition of the two matchings 
corresponding to the two terms. The terms have the same sign if and only if there is an 
even number of loops having one of the two deleted vertices, but not the other, in their 
interiors. 

Definition 1. The Aztec diamond of order n is a planar graph consisting of vertices 

{(2r + 1, 2s) : < r < n, < s < n} U {(2r, 2s + 1) : < r < n, < s < n} (20) 
and of edges 

{{{(2r, 2s + 1), (2r + 1, 2s + 2)}, {(2r + 1, 2s + 2), (2r + 2, 2s + 1)}, 
{(2r + 2, 2s + 1), (2r + 1, 2s)}, {{2r + 1, 2s), (2r, 2s + 1)}} :0<r<n,0<s<n} ^ ' 
in Cartesian coordinates. 



It will soon become apparent that, for, our purposes, the system of coordinates in Figure |l7| 
is more convenient than Cartesian coordinates. It will also become clear why we draw the 
Aztec diamond as if on an infinite square grid tilted 45 degrees from the "natural" direction. 
For now, let us notice that, if we color the Aztec diamond as a checkerboard, all vertices 
on a column have the same color, as do all vertices on a row. Let us also use the system of 
coordinates in Figure |l7| instead of the Cartesian system, and refer to the edge consisting 
of white vertex (xo,yo) and black vertex {x\,yi) as ((#0,2/0)1 (cci , j/i))- 

We can reformulate the rule for comparing terms' signs so that it does not mention loops. 

Lemma 1. Let A and B be two perfect matchings of the Aztec diamond of order n with the 
white vertex at (wq, W\ + d%) and the black vertex at (wq + do, Wi) deleted^. The following 
two conditions are equivalent: 



lr That is, a minor that has all columns of the matrix of which it is a minor, but one, and all rows but 
one. 

2 We shall henceforth use the system of coordinates in Figure Il7] instead of the Cartesian system. 
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one intersection: interior 




Figure 15. Counting intersections 



1. AU B — AD B has an even number of loops containing exactly one of the two deleted 
vertices. 

2. w(A) = w(B) mod 2, where w(T) is the number of edges of the form ((i — l,wi + 
1), (i,wi)), 1 < i < w +d ; ((i,wi + l), (i,wi)), 1 > i < w Q + d ; ((«, u>i (i,wi + 
di — 1), 1 > i < wq, and ((i, wi + di), (i + 1, Wi + d\ — 1)), 1 > i < wq in a perfect 
matching T. 

Proof. The sum w{A) + w{B) is congruent, modulo 2, to the total number of edges in 
AlJ B — A Pi B consisting of a black vertex (x, y) £ D and white vertex (x — 1, y + 1) or 
(x, y + 1). Given a loop £ in AlJ B — Ad B, the number of edges in it consisting of a black 
vertex (x, y) 6 E and white vertex (x — 1, y + 1) or (x, y + 1) is equal to the number of 
times the loop crosses a ray with its end slightly below the deleted black vertex and with 
diagonal direction with respect to the square grid. (See figure |T^.) This number is even if 
the deleted black vertex is in the exterior of the loop, and odd if it is in the interior. The 
same holds for F and the deleted white vertex. Hence a loop has an even number of edges 
consisting of a black vertex (x, y) 6 D and white vertex (x — 1, y + 1) or (x, y + 1) if and 
only if it has exactly one of the two deleted vertices in its interior. We conclude the proof 
by summing over all loops. □ 

It follows that | XX - l) w ^\ — X^tt sg 11 ! 71 ") n™=i a i.ir(i) I) where the sum on the left is over all 
perfect matchings T of the Aztec diamond with white vertex at (wq, w±+di) and black vertex 
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white white white white 

Figure 16. Ordering 



(1,1) (2,1) (3,1) (4,1) 




(1,5) (2,5) (3,5) (4,5) 

Figure 17. Coordinates 



at (wq + do,wi), and |ai,j|™ is the first minor of the Kasteleyn matrix K(G) of the (intact) 
Aztec diamond G where the row corresponding to the white vertex at (wq, W\ + d\) and the 
column corresponding to the black vertex at (wo+do, 101). Thus, if we compute | XX - 1)"^ T '| 
(and this is essentially an enumerative problem, which we shall solve enumeratively) , we 
will know the absolute value of sgn(7r) Yii=i a i,ir(i), which will give us the absolute value 
of the entry of (K{G)~ 1 ) T to be computed, but not its sign. We will find the sign now. 

It is better, for this particular goal, to express the entry of the inverse Kasteleyn matrix 
of the Aztec diamond, not as the ratio of the determinant of the minor of the Kasteleyn 
matrix resulting from the deletion of column i and row jo, divided by the determinant of 
the Kasteleyn matrix, multiplied by (— 1)*°+ J '°, but rather as the ratio of the determinant 
l&ijli to the determinant of the Kasteleyn matrix K(G) — \cn,j\i, where &, j = for 
j ^ jo, h,j = for i ^ i , bi ,j = 1, bij = a,j for i ^ i , j ^ jo- We can then ask 
whether, for a permutation tt of {1, 2, • • ■ ,n — 1}, sgn(7r) ■ Yii ^i^U) nas the same sign as 
(— 1) W ( T \ where T is the perfect matching corresponding to tt. (The answer will be the 
same for all permutations, so we have to ask it for only one permutation (or matching) we 
choose.) If the signs are the same, then E(~ 1 ) ,i ' (T) = (((-^( G )) _1 ) T )io Jo \ K i G )\il if the 
sign are different, then E(-1)"' (T) = -(((^(G))- 1 ) T ) ioJO \K{G)\ n x . 

For our search for the sign, we need to fix an ordering of the black and white vertices of the 
Aztec diamond. The ordering in Figure ^ will prove itself convenient. 

We want to prove that, for any (wo, w\,do,di), 

n 

7T i — 1 



The most straightforward method of proof consists of comparing the signs of J2(— 1)'"'( T ) and 
sgn(7r) n"=i a i,7r(i) for some tiling T and its corresponding permutation tt. Unfortunately, 
this is also quite cumbersome, as which tilings are possible depends on the order of wq, 
wo + do, w\ and w\ + d\. In order not to bore the reader with twenty-four different cases, 
we will proceed explicitly only for the six cases corresponding to do,d\ > 0, which can be 
treated as four different cases. 



15 



2.4.1. Case 1: w\ < wq < W\ + d\ . We choose to examine the matching having /-edges at 
{(i,wi), (i,wi)), wi <i < w Q + d , ((wi,j), (wi, j)), Wi < j < Wi+di, and ((i,wi+di), (i + 
1, Wi+d\ — 1)), u>i < i < wo, where ((xq, yo), {xi,yi)) denotes the edge consiting of a white 
vertex at (xo,yo) an d a black vertex at (x\,y\). All other vertices are covered by \-edges: 
({x,y), (x + 1, y)) for x > y, given that at most one of the two conditions w\ < x < wq + do, 
y = w\ is fulfilled, and ((x,y), (x,y — 1)) for x < y, given that at most one of x = wi, 
w\ < y < w\ + d\ holds, as well as at most one of w\ < x < Wo, y = w\ + d\. (Note that 
we denote by /-edges and \-edges what we called "vertical edges" and "horizontal edges" 
before. This is just a change in notation with the purpose that the fact that we now draw 
the Aztec diamond tilted 45 degrees will not confuse the reader.) It is easy to verify that 
this is a perfect matching. For this matching, (— 1)™( T ) is (— i)« , o+ t "i. 

Each matching corresponds to a permutation of {1, 2, • • • , n}. (We have permutations of 
{1, 2, • • • , n}, and not of {1, 2, • • • ,n — 1}, because we are working with \bij\ 1 , and not with 
a first minor of K{G) = \aij\". The convenience of this choice, which we made without 
justification, will be clear by the end of this paragraph.) The permutation it corresponding 
to the matching we are now examining consists of one cycle. The length of this cycle is 
equal to the number of vertical lozenges, plus one. (This is not true in general. However, 
it is true in cases similar to the one we are currently examining, in which white vertices 
{/l, /2, ■ • • fm} and black vertices {51 }, which would be paired in the form (/», 

1 < i < m in the all-\ matching, are paired in the form (/i,3i+i), 1 < i < m, and vertices 
f m and gi are missing.) Hence sgn(7r) is equal to (—1) to the power of the number of vertical 
lozenges, that is, 2(wo — Wi) + do + d\ — 1 lozenges. 

We now need only examine YliLi Among the edges in our perfect matching, only 

((i, wi + d\), (i + 1, w\ + di — 1)), w\ < i < Wo, correspond to entries equal to (—1) in the 
determinant \b itj \. Therefore UiLi b iMi) = (-l) w °- Wl . 

We conclude that our matching T, for which (-1)">( T ) = (_i)™o+™i-i ; contr ibutcs (_l)2(u>o-u>i)+do+di-i. 
(— l) 1 " 0-1 " 1 , to the determinant. Therefore an arbitrary matching T contributes 

^_^w(T)-(w +wi — 1) _ ^-_-|^2(u>o-"ioi)+<io+<ii-l _ ^_^w -wi {22) 

that is, 

(_l)™(T)(_ 1 j*+di+i ) (23) 

to the determinant. Therefore 

n 

^(-ij-m = ( _i)do + d 1+ i .^ sgn(7r) jj^ 

TV i — 1 

for case 1. 



2.4.2. Case 2: wq > w\ + d\. We choose to examine the matching having //-edges at 
((i,w\),{i,wi)), wi < i < w + d , at ((w + {w , j + 1)), wi + d\ < j < w , at 
{{wi,j),(wi,j)), wi < j < w , and at ((i + 1, w ), (i, w + 1)), Wi < i < w . All other 
vertices are covered by \-edges. The (— 1) W ( T ) of this matching is (— l) dl+1 . The sign of 
the permutation is 

(^Y w o+do-wi-l) + (wo-(wi+di) + l) + (w -w 1 +l) + (w -w 1 ) (24) 
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that is, (-l)<*o+di+i. T h c product IIi=i is (-l) < - w «- ( - Wl + d ^+ 1 ) +< - w °- Wl \ that is, 

(— l) do+dl . Hence each matching T contributes 

(_ 1 )»(T)(_ 1 )do+<i 1 + l (25) 

to the determinant. Therefore 

n 

7T 1=1 

for case 2. 

2.4.3. Case 3: Wq < W\ < w + d . We choose to examine the matching having //-edges 
at ((w , j), (w , j)), w < j < wi + di, at {(i,w ),(i,w )), w < i < w + d , and at 
((wo + do, j + 1), (u> + do + Wo < J < u>i. All other vertices are covered by \-edges. 
The (-1) W ( T ) of this matching is (_i)«>o+tui+i_ Tlle sign f the permutation is (-l) do+dl . 
The product IIi=i is (— 1) M,1 ~ M, °. Therefore each matching T contributes 

(_l)t»(T) . ^_ 1 ^ +di+i (26) 

to the determinant. Therefore 

£ ( -i)«co = .^ sg n W n^ (i ) 

7T 1=1 

for case 3. 

2.4.4. Case 4-' »i > w + d . We choose to examine the matching having //-edges at 
{(w Q ,j),(wo,j)), w <j< w-l +di, at ((i,w ),(i,w )), w a < i < wi, at ((wi, j + 1), (wi + 
1, j)), ui < j < tui, and at ((i, wi + 1), (i + 1, wi)), w + d < i < w\. All other vertices 
are covered by \-edges. Thc (— 1)'"( T ) of this matching is (— l) do+1 . The sign of the per- 
mutation is (-l) d o+di+i p rocmc t f]™ 1 is (-l) do+1 . Therefore each matching T 
contributes 

(_- L yiu(T)(_ 1 ^o+(i 1 + l (27) 

to the determinant. Therefore 

n 

7T i — 1 

for case 4. 



3. An enumerative problem 



Consider an Aztec diamond of side, or order, n, with black vertex (wq + do,wi) and white 
vertex (wo, w\ + d\) missing. Our task in this section is to compute 

B-ir (T) (28) 

T 

where T ranges over all perfect matchings of this Aztec diamond with two missing vertices. 
Enumeratively speaking, we must count the number of perfect matchings, counting as "neg- 
ative" any matching T for which (— 1) W ( T ) = (—1). The function w(T), as defined in the 
previous section, gives, for a perfect matching T, the number of edges consisting of a black 
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vertex [x, y) G D and white vertex (x — l,y + 1) or (a;, y + 1). I? is the union of two set of 
black vertices: E = {(i, Wi) : i < wq + e?o} an d -F = {(«, wi + di) : z < wo}. 

Throughout this section, we will assume do > 0, d\ > 0. At the end it will become clear 
that this assumption involves no loss of generality. One of our tools for attacking the special 
case do > 0, d\ > will be the EKLP Lemma, a classical result in which we will soon 
explain. First of all, we must define two kinds of subsets of the Aztec diamond. They will 
be our intermediate objects of study. 

Definition 2. An n x m black-edged Aztec rectangle with dents at 1 < x\ < x 2 < ... < 
Xm < n+1 is a graph consisting of white vertices 1 < i < n, 1 < j < m, black vertices 

(hj), 1 < i < ti + 1, 1 < j < m — 1, and black vertices (i, m) such that VI < k < m we have 
i =/= Xk- An n x m white-edged Aztec rectangle with teeth at 1 < y\ < y% < ... < y m < n 
is a graph consisting of white vertices 1 < i < 1 < j < m 1 black vertices 

1 < i < n+ 1, 1 < j < m, and white vertices (i, m + 1) s.t. 31 < fc < m with i = Xk- 

Let the number of perfect matchings of an n x m black-edged Aztec rectangle with dents 
at 1 < x\ < x-2 < ... < x m < n + 1 be called D n rn {x\,X2, ...,x m ). What is the number of 
perfect matchings of an n x m white-edged Aztec rectangle Y with teeth at 1 < y\ < y% < 
■ ■■ < Vm < n? Every white vertex (or "tooth") (yi,m + 1) must be covered by a horizontal 
or vertical edge, which will also cover black vertex (y,m) or (y + l,m), respectively. Thus, 
any matching of Y is composed of 

• a matching of a black-edged n x m Aztec rectangle with dents at 1 < x\ < x 2 < ... < 
x m < n + 1) where Xi — yi is either or 1 ; 

• the unique matching of the remaining region. 

Therefore the number of matchings of Y is the m-fold sum 

E n , m (yi,—,y m )= X] D n . m (xi,...,x m ) (29) 

^i—Vi or yi + 1 

where we assume that D n ^ m (x\, x 2 , x m ) — when any two Xj's are equal. 

In the same way that, given D n m , we have found E n ^ m , given E n ^ rn , we can find -Dn. m +i- 
We have again an m-fold sum: 



^2 E ntm (y 1 ,...,y m ) (30) 

Xi<yi<x i+1 

where we assume that E njm (xi, x 2l x m ) = when any two Xi's are equal. 

We will now be able to prove the following by induction. The proof in |3) does not use 
Lemma |^ explicitly. 

Lemma 2 (EKLP Lemma). The number of matchings of an n x m black-edged Aztec rec- 
tangle with dents 1 < X\ < x 2 < ■ ■ ■ < x m < n + 1 is 
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i(m-l) 



D>n,m (x\ 7 X2 ) ■ ■ ■ 3 ^m) 



(m - 1)!! 



where k\\ = 1! 2! 3! ■ • ■ (fc — 1)! fc!. The number of matchings of an n x m white- 
rectangle with dents at 1 < y\ < y 2 < ■ ■ ■ < Vm <n is 



(31) 
Aztec 



E n , m (yi,y2, -,y m ) 



2^ 
(m - 1)!! 



J- 1 



(32) 



For the proof of this lemma, we need a special case of another lemma that we will later 
state in its full generality. 

Lemma 3 (Lemma A (special case)). Let A be an operator carrying polynomials to poly- 
nomials of the same or lesser degree: 



a ofi x u 

a,2.ix 1 + a 2 ,ix 1 + a 2 fiX° 



(33) 



Then 



whe 



A(x j 1 )(x i )\^ = a ,o 



ai,i • • ■ a m _i, m _i 



A(x j )(xi) = a 3tj x{ + a 3 j^ixl 1 + • • • + ajflXi- 



(34) 
(35) 



Proof. The determinant on the left side of ( j34| ) can be obtained from the determinant on 
the right side by elementary column operations. □ 

Proof of EKLP Lemma. The base case (m = 1) is trivial. The inductive step D n ^ m — > E n ^ m 
follows directly from (|30| ) and from Lemma A (special case). 



The inductive step E n 



Dn.m+i requires some more work. 



D n ,m+i(xi, • ■ ■ x m +i) — 



^2 En,m(yi,»>,ym) 
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where B r (y) is the Bernoulli polynomial of degree r, which has the property X^=m 1 



(B q+1 (N) - B q+1 (M))/(q + 1). For a definition of B r (y), consult [Ll|. The only property 



of B r (y) we need to know is that it is a polynomial of degree r whose leading coefficient is 
one. We will adopt, for the sake of convenience, the convention B (y) = y°. We can now 
continue: 



D n ,m+l(xi, ' ' ' Xm+l) — 



_ m(m+l) 


z 


2 










T 


2 



77?!! 



\Bj(x i+ i) - Bj(xi)\™ 

1 B l {x l ) B 2 {x x ) •■• B m ( Xl ) 

1 Bi(x 2 ) B 2 (x 2 ) ••• B m (x 2 ) 

1 B^Xm+i) B 2 (x m+1 ) ■■■ B m (x m +i) 



(37) 



(38) 



where we have pasted a new column onto the left edge of the determinant and a new row 
onto the top edge, and then added the first row to the second one, the second to the third 
one, and so on, in succession 



raw 



by Lemma A (special case). 



J- 1 



m+1 
1 



(39) 
□ 



Let us now count matchings with weights 1 and (—1): a matching may count as one matching 
or as minus one matchings. Consider an n x m black-edged Aztec rectangle with dents at 
\ < x\ < x 2 < ■ ■ -Xm-i < n + 1, and, in addition, a dent at wq. Let us multiply the total 
number of matchings by (—1) if there is an odd number of x^s smaller than wq. (Here we 
are counting either all matchings as positive matchings or all as negative matchings.) What 
is, then, this total, weighted number of matchings? Conveniently, it is 



(m-l)H 



IV, 



i-i 



x\ 
x 2 

3-1 
X m-'i 



(40) 



We have just taken ( |3l| ) and shifted the w 3 row corresponding to the dent at wo as many 
positions upwards as there are x^s smaller than uuq. 

What is the weighted number of matchings of a white-edged n x m Aztec rectangle Y with 
teeth at yi, y 2 , ■ ■■ y m -i and a black hole at (wo, m)l (The weight of each matching is (—1) 
to the power of the number of its edges consist of a black vertex (i, Wx) and white vertex 
(i — 1, w\ + 1) or (i, Wi + 1), where 1 < i < wq + do.) A cursory examination makes clear 
that a matching of Y will have weight 1 if and only if the black-edged nxm Aztec rectangle 
the matching outlines has an even number of dents with indices lower than w. Thus the 
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weighted number of matchings is 



which is the same as 



5^ I in - 1)! 

i=ViOryi + l 



i-i 



i-i 

c m-l 



2^ 
(m- 1)!! 



yr 1 + (i/i + 

yr 1 + (?/2 + ly- 1 



vL-i + (ym-i + ly- 1 



(41) 



(42) 



Here we face a difficulty. We would like to use Lemma 0. Unfortunately, we have w 3 on 
the top row. What we need is to find something which is transformed into u>^ -1 by the 
linear operations taking a row of the form x 3 ^ 1 to x 3 ^ 1 + (x + l)- 7 ' -1 . First of all, we need 
a stronger version of Lemma [}[ 

Lemma 4 (Lemma A, stronger version) . Let A be an operator carrying polynomials to poly- 
nomials of the same or lesser degree: 



+ aifiX 



a.2,2X 2 + a 2 ,ix 1 + a 2 fiX° 



(43) 



Then 



Kk<L 



OO.O—Om-l.m-l 



Kk<L 



(44) 



for any l i; bk i, 1 < k < £j, 1 < i < m, where 

A(x r )(xi) = a jd xl + ajj-ixi' 1 + 



(45) 



Proof. The special case followed from the fact that the same sequence of column operations 
transforms the row A(x 3 ~ 1 )(y) into y 3 and the row A{x 3 ){z) into z 3 , for any values y,z. 
Therefore the same sequence transforms the row a-A(x 3 )(y)+b-A(x 3 )(z) into a-y 3 +b-z 3 . □ 



We need afc's such that J2k a kA{x 3 ){vk) = Wq for < j < m, where A is the operator 
taking x 3 to x 3 + (x + I) 3 . 

Definition 3. A is the operator taking x 3 to (x + l) 3 — x 3 . 
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Then A = (21 + A), where / is the identity operator, and we have, formally, 



(27 + A)-'=i-(/+^)-' 



-i-D-nf <46) 

3=0 



Since A J vanishes on polynomials of degree smaller than j, and we are dealing with poly- 
nomials of degree at most to — 1, we can use just the first m terms of the series: 



m— 1 

(l/2)£(-l)'(A/2)' (47) 

3=0 

The reader can easily check that, for any x k with k < to, 

1 m_1 A 

(2J + A)((- £(-!)'(-)')(**)) (48) 



gives a; fe plus a constant times A m x k , and, since m > k, A m x k is zero. One can also check 
the same for 



m-l a 

«2 £(-l) J '( T ) J ')(* fc ))(2J + A). (49) 



It is quite convenient that (21 + A) and 

m— 1 

v 2 



m-l a 



j=0 



commute, and that their composition is the same as the identity operator for the domain 
we are interested in. We will use the shorthand (21 + A) -1 for 



m—l 

((-]T(-ir(A/2>0(x fe )) (51) 



3=0 

without any compuctions. 
If we define do, a\, • • • , a TO _i by 

rn — 1 m—l 

v 2 



rn — 1 m — 1 

(- £(-iy(A/2)^)(/0r)) - £ */(* + »), (52) 



j=0 i=0 
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then 



nt — L in — j. 

£ ^((2/ + A)(xi)(w +i)) = (((- £(-l)'(A/2)')(2I + A))(^'))( Wo ) 

= (z J )(wo) 



for < j < m, as we desired. 
Hence, by Lemma A, 



= w i 



w, 



j'-i 



yf 1 + (2/1 + 

yf 1 + (2/2 + 1)'"" 1 

l4"-l + (l/m-1 + l)'- 1 



is equal to 



2 m . 



(lEL'o^-i)^!)")^- 1 )^) 

2/i 

j'-i 

2/2 



2/m-l 



or, in shorthand, 



((2/ + A)- 1 (^- 1 ))( U ; ) 



m 

2/2 



Vm-l 



Now, for every polynomial p, 



p( Wo ) = (((/+Ar- i )( P ))(i). 



Hence we can write 



2 m . 



((/ + Ar°- 1 ((2/ + A)- 1 (^'- 1 )))(l) 



j'-i 
2/i 

j-i 

2/2 



Vm-l 
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and eliminate all terms of degree m or higher. Hence the weighted number of matchings 
of a white-edged n x m Aztec rectangle with teeth at 2/1,2/2, • • -J/m-i and a black hole at 
(wo,m) is 



m(m + 1) 



(m-l)H 



(( / + A )»o-i(( 2/ + A)- 1 (^- 1 )))(l) 



7— 1 
2/2 



j'-i 

Vm-l 



(59) 



What would be the weighted number of matchings of a black-edged n x (m + 1) Aztec 
rectangle with dents at x\, X2, • ■ • x m and a black hole at (wo 7 m)? This is the sum we have 
to simplify: 



2^ 
^ (m-l)H 

Xi<Vi<x i+1 



((I + Ar°- 1 ((27 + A)- 1 (^'- 1 )))(l) 



j'-i 
J/i 

j'-i 
2/2 



(60) 



m(m + l) 

This is equal to 2 ( m _ 2 1 v ! times 
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E 

Xi<Vi<x i+1 



((/ + A) W °~\(2I + A)~ 1 (a;- 7_1 )))(1) 
j'-i 

i-i 

2/2 



j'-i 

2/m-l 

((/ + A)»»- 1 p + A)- 1 (^- 1 )))(l) 
a^ -1 + (an + ly'- 1 + --- + (x 2 - l)'- 1 

A^ 1 + + ly- 1 + ■ ■ ■ + (x 3 - 

((/ + A)-o-i((27 + A)- 1 (^- 1 )))(l) 
4(5,(^2) -S^xO) 
^•(a*) -S j ( a;2 )) 

1 jrr-Sj-i^i) 

((/ + A)-«- 1 ((27 + A)- 1 (^" 2 )))(l) 

j^Bj.^) - Bj.^X!) 



jzr[Bj-i(x m ) — Bj_i(x m _i) 

((/ + A)-«- 1 ((2/ + A)- 1 (^'- 2 )))(l) 

1 jhBj-^i) 

Tzr-Sj-i^) — Bj_i(xi) 







((/ + Ar°- 1 ((2/ + A)- 1 (^"- 2 )))(l) 



=1^-1 (*i) 



/— l Bj—i (x m ) 



Since {((/ + A) tu °- 1 ((27 + A)" 1 (a; : '- 2 )))(l)} j= 2,3,... m+ i is a linear combination of 
the form {fc- J_2 }j = 2,3,... m +i, 1 < k < m, it is enough to show how to simplify 



¥~ 2 

1 -pjBj_i(xi) 
1 O^-i (x 2 ) 

1 lh B j-l{Xm) 
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This happens to be easier than one would expect. As a particular case of 




B q+1 (N) - B q+1 (M) 
q + l 



we have 

= kB 3 (k + i)-B 3 (k)). 

Therefore 



y- 2 

1 J^Bj-^X!) 
1 i S,.!^) 



1 j=T B j-l( x m) 



1-1 ^(B^k + ^-B^k)) 
1 -pjBj_i(xi) 
1 j^jBj-i(x 2 ) 

(k + ly- 1 - y- 1 

i-i 
j-i 

2-2 



A(^'- 1 )(fc) 
x\ 

2 



m! 



(68) 



(69) 



(70) 



Thus, we now know that the sequence of elementary column operations we have to ap- 
ply to the matrix having x\ on its lower rows in order to make it into a matrix having 
{1, jBi(xi), \B-2ixi), ■ ■ ■ ^B m (xi)} on its lower rows transforms the row {A(xi~ 1 )(k)}™ +1 
into the row {0, k°, k 1 , ■ ■ ■ k m }. What row would be transformed into the row 

{0,((I + A)'"»- 1 ((2/ + A)- V))(l), ((/ + A)— 1 ((27 + A)- 1 ( 2 ; 1 ))(l), 
((/ + Ar°- 1 ((2I + A)- 1 (x m ^ 1 ))(l)}? 

Let us express (I + A) W °~ 1 ((2I + A)~ 1 )(f)(x) (/ any polynomial of degree at most m — 1) 
in the form Y^k=0 a kf(x + k) for some numbers oi, 0,2, ■ ■ ■ a m . If rows J/i, t/2, • • • 2/r are 
transformed into rows Zi, Z2, • ■ • z r , respectively, then the row Y^i=i biVi must be transformed 
into the row Y"^_i biZ;. Therefore the row 

{Efco 1 a k {A{xi- l )){k + 1)}™ +1 is transformed 

into 

m— 1 m— 1 m— 1 

{0, ^ a fe (fc + 1)°, «fe( fc + i) 1 . ' ' ' . E a *( fc + i)" 1 " 1 }' 

fc=0 fc=0 fe=0 

which is the same as 

{0,(7 + Ar°- 1 ((2/ + A)- V))(l), (7 + A) W °-\(2I + A)~ V))(l), 
(I + Ar°- 1 ((2I + A)- 1 (x m - 1 ))(l)}. 
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Now, what is J^kLo a fc(A(a; J ))(fc + 1) for < j < ml It is 

((/ + Ar- 1 ((2/ + A)- 1 (A(^))))(1). 



Hence the weighted number of matchings of a black-edged n x (m + 1) Aztec rectangle with 
dents at x\, ■ ■ ■ , x m and a black hole at (wq, m) is 



((/ + A) W °- 1 (2I + A)- 1 A)(^'- 1 )(1) 



3-1 

3-1 
C 2 

ri-1 



(71) 



In the same way we have arrived at this result, and using it as a base case for induction, it 
is easy to prove the following. 

Lemma 5. The weighted number of matchings of a white-edged ri x (m+di) Aztec rectangle 
with teeth at y\, ■ ■ ■ y m +di-i an d a black hole at (wq, m) is 



(m + d 1 )(m + d 1 + l) 



(-1) 



di ■ 



(m + di- 1)!! 



((/ + A) W °- 1 (2I + A)-( dl+1 )A dl )(^'- 1 )(l) 



3-1 

j-i 
vi 



3-1 
Vm+di-l 



(72) 



and that the weighted number of matchings of a black-edged n x (m + d\) Aztec rectangle 
with dents at x\, X2, ■ ■ ■ ,X m +d\—i an d a black hole at (uio,m) is 



(-IF- 



( m +d 1 )( m +d 1 -i) 



(m + di - 1) 



((/ + A) W «- 1 (2I + A)- dl A rfl )(^ _1 )(l) 



3-1 
3-1 



3-1 
Vm+di-l 



(73) 



Notice that we have used the fact that (21 + A) -1 and A commute. They do so because 
(2/ + A) -1 , on the domain of polynomials of degree lower than a given bound, is shorthand 
for a finite sum of powers of A. 

We can now attack our main objective, namely, the enumeration of matchings of an Aztec 
diamond with a black hole at (wq + do, w±) and a white hole at (wq, w\ + d\), where some 
matchings are counted as negative matchings. Let us first consider the special case wq = 1. 
What happens at white and black vertices (x, y) with x = 11 The hole at (l,w\ +d±) forces 
a zig-zag pattern covering all those vertices. (See figure |l8|.) None of the edges covering 
these vertices counts towards the weight of the matching. Thus, the weighted number of 
matchings of an Aztec diamond of order n with a black hole at (1 + do, w\) and a white hole 
at (1, iui +d\) is the same as the weighted number of matchings of a (n — 1) x n white-edged 
Aztec rectangle with teeth at {1, 2, ■ • ■ n} and a black hole at (do,w±). By Lemma |J this is 
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equal to 



(-1)' 



2^ 
(n-1)!! 



((/ + A) d °- 1 (27 + A)-("- ,Ul + 1 )A™- M,1 )(^- 1 )(l) 
(n - ly- 1 



(74) 



The top row is a linear combination of rows of the form {A k (x j - 1 ){l)} 1 1 < j < n. All 
such rows with k > (n — 1) vanish. All such rows with k < (n — 1) are linear combinations 
of the rows 

{f 1 }, 1 < j < n, 1 < i < n - 1, 

Since an n by n determinant whose last n — 1 rows are {2 J }, • • • {(ri — and whose 
first row is one of ^ vanishes, we can discard all terms of the form 

{A fc (^- 1 )(l)},l<j<n,A : <(n-l) 



from the top row of the determinant in ( |74| ) . Hence we must consider only the term of the 
form 

in the top row of the determinant. Thus ( [74] ) is equal to 

nCn+i) 

(-1)™ -r^yyy • (C ■ ■ (n - 1)!!) 

= C?-(-l) ,i,I + 1 2 M ^. 

times the coefficient of A™ -1 in 

(I + A) d °- 1 (2I + A)- (n+1 - Wl ) A n ~ Wl . 

Since, in the task of finding this coefficient, A plays a purely symbolic role, we might as 
well use a symbol which does not denote an operator, such as x. It will also be convenient 
to use the following notation, due to Richard Stanley. 

Definition 4. The coefficient of in the formal power series f(x) on x is denoted by 
[x j ](f(x)). 



2,x 



Our task is to compute [ar n_1 ]((l + x) d °- 1 (2 + x )-(n+i-wx) X n- Wl _ We have 

[x n - 1 ]((l + x) d °- 1 {2 + x )-(.n+l-wx) x n-w^ _ + a;)*>-l(2 + 

= [(2^r- i ]((2+ 2y r- i ii±|g^ ) 

= a-^ 1 -^ ■ [y u ' 1 - 1 ]((2 + 2yr- l (1 +^ ) °~ 1 ) 

= 2-< Wl - 1 ) ■ 2— 1 2-« • [^-^((l + ?/ri -i ( 1 +^°- 1 

= 2-" • br^Ka + ^-i ^+y)'"" 1 ) 

Now, for any formal power series f(y) on y, 

Wl —1 



W1- -L / _ 1 \ 

^-^((i+vr- 1 •/(»))= E \ ji^k/go) 

-i— n \ J / 



3=0 

= E 1 ( K " 1 : j)+i )M(/(.)) 

i=o ^ > 

wi-i j (75) 

= E ((^'"^ n _ z . )j+1 )) • (WKfiv)))) 

3=0 [ ' 

= S «[»" 1 K (i :^ ) j-n » ■ di^KS ■ /(»)))) 



]((t^-)-/(t^-)) 
1 — z 1 — z 



Therefore 

( CI i x) d ^ 1 (2 + x)-^ 1 -^) = 2- n ■ [y^Kil + ^- 1 il±M^ )( by (D) 

_2z_\d — 1 
-) 



= 2~ w • 




-\K(1 


= 2"™ • 


[z W1 ' 


%. 


= 2"" • 




%. 


= 2"" • 


[z w ^ 


%. 


= 2"" • 


[z Wl ' 


" x ]((l 



_5_ W 



1-Z (1 + ^ 

1 2z _ (76) 



1 -z' 

(^j_£)*-i .(i-z)™) 



1 ,l + ^d -l 



2!) 



An expression of this form is called a Krawtchouk polynomial (||, p. 130). 
We have just proven 

Lemma 6. The weighted number of matchings of an Aztec diamond of order n with a black 
hole at (do + and a white hole at (l,w± +d±), do,d\ > 1, is 

. ^"^((l + zf"- 1 ■ (1 - z )(n-l)-(*-l)) . <2r^. 



We can now work on the general case. Every matching of an Aztec diamond of order n with 
a black hole at (wq + do, w±) and a white hole at (wo, w\ + e?i) can be subdivided into one, 
and only one, pair of matchings of the following form. 



1. The first item of the pair is a matching of a white-edged n x (w\ + d% — 1) Aztec 
rectangle with dents at yi, j/2> ' " ' Dwx+dx-i and with a black hole at (wq + do, u>i). 

2. The second item of the pair is a matching of a white-edged n x (n — (w\ + d\ — 1)) 
Aztec rectangle with dents at Z\,z%, ■■ ■ £ n — (iui+di— !)• 



The numbers yi, y 2 , ■ ■ ■ ^+4-2, Z\, Z2, ■ ■ ■ z n -( W \+di-i) an d w o are all distinct and cover 
all of the interval {1, 2, • ■ • , n}. Thus every matching can be described as a partition of 
{1,2,--- , n} — {wo} into two subsets, a matching of an Aztec rectangle with the first 
subset as its set of dents, and a matching of an Aztec rectangle with the second subset 
as its set of dents. As we stated in the previous section, we will weigh each matching 
T by a factor of (-1) W< - T \ where w(T) is the sum of the number of edges of the form 
((i — 1,W\ + 1), (£, Wi)) or {(i,wi + l),(i,wi)) for 1 < i < Wq + do and the number of 
edges of the form ({i,w\ + d\), (i,w\ + d\ — 1)) or ({i,w\ + d\), (i + l,wi + d\ — 1)) for 
1 < i < wo- Thus the weight w(T) of a matching is equal to the weight of the matching 
of the n x (wi + d± — 1) white-edged Aztec rectangle which it induces, plus the number of 
teeth of this Aztec rectangle whose indices are lower than wq. 



Hence the weighted number of matchings of an Aztec diamond of order n with a black hole 
at (wo + do, w±) and a white hole at (wq, w± + d±) is equal to the sum of 



(^ 1 +d 1 -l)(u 11 +d 1 ) 



(tflx + di -2)!! 
((7 + A)™"^- 1 ^! + A)-*A*- 1 (3;j- 1 )))(l) 



y{ 

i-i 



i-i 

Vw-i+dx-2 

, (■n.-(-w 1 +d 1 -l))((n-(u, 1 +d 1 -l)) + l) 



■/J- 1 



i— {wx+di— X) 



(77) 



(n- (wi + di - 1) - 1)!! 
over all partitions of {1, 2, • ■ • , n} — {wo} into two sets 

{Z/l> V2, ■ ■ ■ y wl +d!-2}, { z l, z 2, ■ ■ ■ £ri-(ttJi+di-l)}> 

where y x < y 2 < ■ ■ ■ < y Wl +d x -2 and z\ < z 2 < ■ ■ ■ < 2r„_( Wl+c2l _ 1)) and t(w ,yi, 
equal to how many of t/i, y%, ■ ■ ■ yu are less than Wq. 



,Vk) is 



:!() 



We may dispose of the inconvenient t(wo, yi, ■ ■ • , ywi+di-2) by expressing the same result 
as the sum of 

((/ + A) ro °+ d °- 1 ((2/ + A)-*A*- x (a^- a )))(l) 



(-1)** 



(wi + di-2)!! 



u v wi+dl -i,w v w 1 +d 1 -l 



2" 



(r.-( mi +d 1 -l))(( re -( mi +d 1 -!)) + !) 



(n - (ti7 X + d x - 1) - 1)!! 



(78) 

over all partitions of {1,2,- •• , n} into two sets {ui,-- - ,v Wl +d!-i}, Z\, ■ ■ ■ z n -i W i+dx-X)i 
where V\ < V2 < ■ ■ • < v Wl+ d 1 -i and z\ < z 2 < ■ ■ ■ < z n -(w 1 +d 1 -i), and where Sij is equal 
to 1 for i = j, for i ^ j. 

The following lemma follows immediately from Laplace's development of a determinant 
(||, pp. 22-25). 

Lemma 7. Let us have a determinant \bi,j\™ 1+m2 > where bij — Cij for j < m 1; h t j — 
{-\) % - x dij- mi for j >m 1 . Then 

\bij\T 1+m2 = (-l)L^-^J • J2 • \d*ui\T (79) 

where the sum is over all partitions of {1, 2, ■ • ■ , mi + TO2} into two sets {xi, X2, ■ • ■ , x mi }, 
{yi, 2/2, ■ • ■ , Vm 2 }, where x ± < x 2 < ■ ■ ■ < x mi and y± < y 2 < ■ ■ ■ < y m2 , and where [z\ is 
the largest integer less than or equal to z. 

We can now express our sum over partitions as a determinant. The weighted number of 
matchings of an Aztec diamond with two holes is 



( roi+ d 1 -l)( mi +d 1 ) („- (mi+ d 1 -l))((r.-( mi +d 1 -!)) + !) 

(-1)^- . - f_i)L^ i J-L^4 £l J \ d - (80) 

[ ' {w 1 + d 1 -2)\\ ( n -{w 1 + d 1 -l)-l)\\ [ ' ' [ ' 

where 

• di,i = for all 1 < i < n + 1, i ^ Wo + 1, 

• d Wg+ i.i = 1, 

• d ld = ((/ + A) Wo+d °- 1 ((2I + A)- dl A dl - 1 (a; J '- 2 )))(l) for all Kj< w 1 +dt, 

• di j = for all w± + di < j < n + 1, 

• dij = (i - l) j ~ 2 for i > 1, 1 < j < wi + d\, 

• d- tj = (-l)*" 1 ^ ~ iy-(^+ d ^+ 1 ) for i > 1, j > wi + di 

The task ahead is to compute the determinant |djj|? + . For convenience, we will refer to 
it as D(wo, do, w\, d\). From (|8fj| ) and Lemma it follows that 

(™ 1 + d 1 -l)(v, 1 +d 1 ) („-( mi +d 1 -l))((n-( TO1 +d 1 -l)) + l) 

ua,*.-.,*) - «-!>'■ (Bi - - m in _ iwi - 1 _ T) _ 1)n (-Dtfj-t^j)- 

(-1)™ 1+1 • [^ 1_1 ]((1 + zf ' 1 • (1 - z )(»-i)-(do-i)) . 2 Ili V Li . 
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We will reduce the general case to the special case u> = by expressing D(wq, do, Wi, d\) — 
D(wo — 1, do, wi, di), as the product of D(l, do + w — 1, wi, d\) times something else. 

If we take the determinant |dj ; j|™ +1 for D(w — l,do,Wi,d\), and add one to the bases of 
all powers, we obtain 

D(w -l,d ,w 1 ,d 1 ) = \g iJ \^ +1 (81) 

where 

• gi y i = for all 1 < i < n + 1, i ^ w ; 

• g wo ,i = i; 

• 5ij = ((/ + A) M,0 +' i( '- 1 ((2/ + A)- dl A dl - 1 (a;J- 2 )))(l) for all 1< j < wi + di ; notice 
how we have raised the exponent of (I + A) from Wo + d n — 2 to w + d — 1; 

• gij = for all w\ + di < j < n + 1; 

• gjj = z J ~ 2 for i > 1, 1 < j < wi + di, 

• 5ij = (-l)«- 1 iJ"-(^i+ d i+ 1 ) for i > 1, j > Wl + di, 

The n-tuple |5n+i,jl2 +1 i s a linear combination of the n-tuples |<7i,j I2 +1 ' 2 < i < n + 1 and 
of the n-tuple {l}2 + \ which is what gij would be if the pattern for {g n ,j}2 +1 1 {dn-~L,j}2 + > 
• • ■ {.92j}2 +1 were continued. 

Lemma 8. Let a k be the coefficient of x k in ((x - l) Wl+dl - 1 (x + 1 )™-(™i+<ii)+i). Then 

n 

^a k {k + iy- 2 = for 1< j <wi+di, (82) 

fc=0 

n 

a k {-l) k {k + iy-( w i+ d i +1 1 = for w 1 +d 1 + l<j<n+ 1, (83) 

fc=0 

Proof. 

A Wl+dl - 1 x : >- 2 = for 1< j < wi + di (84) 

implies 

A ™i+di-i( A + 2 /)n-(«'i+di)+i a j-2 = for 1< j < wi + di (85) 
If we take the value at x = 1, we obtain 

n 

J2a k (k + iy- 2 = for 1 < j < wi+di. (86) 

fc=0 

Similarly, 

A „-( wl+dl ) + i xj - {wi+dl+ i) = for wi + di + 1 < j < n + 1 (87) 

implies 

(A + 2I) Wl+dl -\-K) n - {wi+d ^ +1 x 1 - {wi+dl+1) = for wi + di + 1 < j < n + 1. (88) 
If we let -ff be the operator taking to (x + l)- 7 ', we can write 

(7 + tf^i+^-^J - fl-jn-Ctoi+dO+igj-^i+di+i) = for wi + di + 1 < j < n + 1 (89) 



:s2 



Clearly the coefficient of H k in 

f^j _|_ jjyui+di— 1 _ jjyi-(wi+di)+l 

is equal to (— l) fe times the coefficient of H k in 

_ jjyvi+di-l _|_ + l 

which is equal to (— i^i+di- 1 times the coefficient of H k in 
that is, afe. Hence 

n n 

(.ijtBi+di-i ^(_i)* afc (jfe + = ^((_i)^i+rfi-i(_i)fe afe )( fc + ^j-^i+di+i) 

= _|_ jj-^wx+dx-l _ jjyn-(wi+di)+l x j-(wi+di+l) ^ 

= (A + 2I) Wl+dl ^ 1 {— /\) n ^^ Wl+dl ^ +1 x^^ Wl+dl+ ^ 
— (A + 2/)"' 1+ ' il ~ 1 A"~("' 1+cil - )+1 2; : '~( M,1+ ' il+:L ' ) 
= for wx + dx + 1 < j < n + 1. 

Therefore 

^ afe(-l) fe (A: + l)J-(«>i+<ii+i) = for wi + di + 1 < j < n + 1. 

fc=0 

□ 

Hence, if we add a n times the nth row, a n -\ times the (n — l)th row, . . . 02 times the 
second row to the bottom row of , we obtain the row where 

• /i = [^"-^((a; - l) Wl+dl - 1 (a; + 

• /ij = (-1) • for 1< j < wi + di, 

• ftj = (-1) • (-lj^i+di-i fo r j > Wl + d x . 

After switching signs and shifting the bottom row of D(wq — 1, do, w\, di) (now {/ij}™ -1-1 ) 
to the second-to-topmost place, we obtain 

D(«, - l.do.wi.di) = (-1) ' (-l)^*" 1 • • |fc„^ +1 , (90) 

where 

• h,i = for all 1 < i < n + 1, i ^ Wo + 1, « ^ 2; 

• fe2,i = (-1) • (-l) 1 "^- 1 • [a^ " 1 ]^ - l)"^ 1 " 1 ^ + i)™+i-(™i+<ii)) ; 

• feu>o+i,i = i; 

• k hj = {(I + A) w "+ d «- 1 ((2I + A)- dl A dl - 1 (^'- 2 )))(1) for all 1< j < w x + di; 

• fcij = for u>i + di + 1 < j < n + 1; 

• fejj = (i - l)^ 2 for i > 1, 1 < j < wi + di; 

• fcij = (-l) 4 (i - for i > 1, j > Wl + di + 1. 

We multiply the n — (idi + di) + 1 rightmost columns by (— 1), obtaining 

Z)(tuo - l,d ,wi,di) = , (91) 

where 
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h,i = for all 1 < i < n + 1, i ^ wq + 1, % ^ 2; 

fe'.i = (-1) • (-l)^i +^1-! • [^"^((a; - l) ,Ul + <il - 1 (a; + i)"+i-(«>i+<ii)) ; 

^u.o+1,1 = i; 

= ((/ + A) u *>+*- 1 ((2I + A)- dl A dl - 1 (^'- 2 )))(l) for all 1< j <wi+di; 
Zij = for wi + d\ + 1 < j < n + 1; 
Zi,j = (i - 1) J ~ 2 for i > 1, 1 < j < wi + di; 
l-j = - lp'-^i+di+i) for i>l,j>w 1 +d 1 + l. 



Therefore 



D(w ,d ,wi,di) - D(w - l,do,wi,di) = |ri,j|" +1 , (92) 



where 



• r,,i = for all 1 < % < n + 1, i ^ 2; 

• r 2 ,i = (-l)«'i+ d i- 1 • [a;"' - 1 ] ((a; - l) Wl+dl - 1 (a; + i)™+i-(«>i+<ii)) ; 

• n'j = ((/ + A)"'° +<i °- 1 ((2J + A)- dl A dl - 1 (a;J- 2 )))(l) for all 1< j < w x + d i; 

• rij = for !«! + di + 1 < j < n + 1; 

• r it j = (i - l) j ~ 2 for i > 1, 1 < j < w\ + d\\ 

• r Lj = (-ly-^i - iy-(^+^+ 1 ) for i > 1, j > Wl + d 1 + 1. 

This is equal to 

r 2 ,i-D(l,itfo + do-l,u>i,di) = (_ip+<*i-i . [a^ " 1 ]^ - l) Wl+dl -\x + i^+M^+di)) 

D(l,w + d - l,ioi, di) 

= ^"^((l - aj)™i+<*i-i(i + ^("-(^i+di-i))) . 
D(l,w + d - l,ioi,di) 



(93) 



Therefore 



"Jo — 1 

D(w ,do,wi,di) = ( ^2 D U + l,d ,wi,di) - D(j,d ,ioi,di)) 
+ D(l,do,toi,di) 

wo — 1 

= ]T ([^](1 - ^ + ^""(-i+^-l)) . £)(1, j + do,U>l,di) 

+ D(l,do,tui,di) 

t up — 1 

= (t^'K 1 " ^r i+dl_1 (l + a;)"-^ 1 ^ 1 " 1 )) • D(l,j + d ,w 1 ,d 1 ) 

3=0 

(94) 
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By (|S0|), it follows that the weighted number of matchings XX ~ l)™^ of an Aztec diamond 
of order n with a black hole at (wq + do, it?i) and a white hole at (wq, w\ + d\) is equal to 

(u, 1 + d 1 -l)( W1 +d 1 ) (n-(v, 1 +d 1 -l))((n-(u 11 +d 1 -!)) + !) 

1 ' (w 1 +d 1 -2)\\ (n-(ioi + di-l)-l)!! v ; 

w -l 

£ - + x)"-^^ 1 - 1 )) • D(l,j + do,w u d x ) 

■ wq- 1 

= (([^](i-^) iui+di_1 (i+^) n_(iui+di_1) )- 

.7=0 

( „_ (mi+dl _i ))(( „_ (mi+dl _i )) + i ) 

vv y (wi+di-2)!! (n - (mi + di - 1) - 1)!! 
(-ljL^J-L^J.Dfl.i + do.iDi.di)) 

The term within parentheses including D(l, j + do, w±, di) is equal to the number of match- 
ings of an Aztec diamond of order n with a black hole at (J + do + 1, ioi) and a white hole 
at (1, wi + di). By Lemma this number is equal to 

(_l)^i+i . 2 Ili2 ^ 11 ■ [s u,1-1 ]((l + zy +d °- 1 ■ (1 - 
Hence XX - l) to( ' T ' ) is equal to 

t up— 1 

• 2-" ■ J2 (i xj W - x) Wl+dl - 1 (l + a;)»-(wi-HJi-i)) . 

j=0 

[^ 1_1 ]((1 + z)^"- 1 ■ (1 - 2 )(»- 1 )-(i+ d "- 1 ))) i 

From this and from (|l|) the result we have sought follows immediately. 

Proposition 9. The entry in the inverse of the Kasteleyn matrix of an Aztec diamond of 
order n corresponding to a black square at (vjQ + d ,wi) and a white square at (wo,wi+di), 
d ,di > 0, is 

Wq— 1 

rydo+dx+m . 2 -n . \^ (M'l((l - x) m+dl - 1 (l + x) n -( wl+dl -V) ■ 

U ( 95 ) 

where [x J ](p(x)) is the coefficient of x 3 in the polynomial p(x) . (Alternatively, this can be 
called the value of the coupling function of the Aztec diamond of order n at the black hole 
(wq + do,Wi) and the white hole (wo,w\ +d\).) 

In order to deal with the cases do < and di < 0, we merely need to flip the Aztec diamond 
so as to make do,di > and compute the weighted number of tilings in the manner we 
have described. Of course, we have to account for the fact that the weighting has to be 
computed differently. For do < 0, we also have to express D(wo, do, Wi, di) as 

n 

- ^2 D U + l,d ,u>i,di) - D(j,d ,wi,di) 
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and not as 



Wq— 1 

D(l,do,wi,di) + D(j + l,d Q ,wi,di) - D(j,d ,wi,di)) 
i=i 

as we did in (p4|). These are the only two details worth mention in the otherwise trivial 
derivation of the following result from Proposition |9|. 

Corollary 10. The coupling function of the Aztec diamond of order n at the black square 
(wo + do, wi) and the white square (wo,wi + di) is 

w -l 

U (96) 

[« W1_1 ]((1 + z) 3+d °- X ■ (1 - z )("-l)-(i+<io-l)-j) 

for do > 0, 

n 

2 -» . (_ ^ ([^']((1 _ x )»i+di-l(l + a; )n-(»i+di-l)) . 

[^"^((l + z )i+<*o-l . (1 _ z )(«-l)-(i+d -l))-j) 

/or d < 0. 

When we take a minor of the inverse Kasteleyn matrix, the factors (— l)«o+«i+«>i j multiplied, 
give the same product in every term of the expression of the minor in a form such as 
sgn(7r) Y\i=i c i.-n(i)- Thus we can leave them out, and our main result follows. 

Theorem 11. The probability of a pattern covering white squares v%,V2,- ■ and black 
squares Wi,W2, ■ ■ ■ Wk of an Aztec diamond of order n is equal to the absolute value of 

l c («i>^-)li, J -=i,2,...jfc- 
The coupling function c(v, w) at white square v and black square w is 

2~ n kr(j, n, Vl - 1) kr^ - l,n - 1, n - (j + xh - x t )) 

3=0 

for xii > Xi and 

n 

-2~ n Y kr 0'= n ' Vi ~ l ) kr (y'i - 1, n - 1, n - (j + xii - Xi)) 

j=Xi 

for xU < Xi, where {xi,y{) and are the coordinates of v and w, respectively, in the 

coordinate system in figure [T^J, and the Krawtchouk polynomial kr(a, b, c) is the coefficient 
ofx a in (l-x) c -(l + x) b ~ c . 




Figure 19. Inverse Kasteleyn matrix entries as a function of wq, w% for 
d = 1, di = 2, n = 40 



4. In perspective 



Proposition |9| is valuable in itself, in that it gives us an efficient algorithm for computing 
an arbitrary entry in the inverse Kasteleyn matrix of the Aztec diamond. Figures [l^ to 
p4| show the absolute value of the entry as a function of wq and Wi , for fixed do and di . 
As we showed in section 2, given a pattern consisting of k vertices, we can compute the 
probability of its occurence at any point in a Aztec diamond of given order by computing 
(fc/2) 2 entries of the inverse Kasteleyn matrix. Thus, for a fixed pattern, the time required 
for computing its probability is equal to a constant times the time required for computing 
an entry of the inverse Kasteleyn matrix using Proposition ^. Whether computation time 
grows quadratically on the order of the Aztec diamond, or somewhat faster, depends on 
whether multiplying integers is assumed to take constant time. What is clear is that we 
now have an algorithm that is much more efficient than computing the entries of an inverse 
Kasteleyn matrix by actually inverting the matrix or computing minors. 

Proposition ^| gives us an expression that is more closed than an entry in the inverse of a 
Kasteleyn matrix. What do we mean by this? There are few tools available that would 
allow us to obtain asymptotic expressions for a sequence of entries in a sequence of inverses 
of arbitrary matrices. For finding the asymptotics of sums such as (^5|), however, there 
are many well-developed analytical techniques. At the time of this writing, Henry Colin 
is working on some minor problems involved in applying the saddle-point technique to the 
asymptotics of (^) . Once he superates these difficulties (something that seems to be about 
to happen), the goals set in the introduction will have been achieved completely. 
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Figure 20. Inverse Kasteleyn matrix entries as a function of wq, w\ for 
d Q = 1, di = 2, n = 80 
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Figure 21. Inverse Kasteleyn matrix entries as a function of wq, w± for 
do = 1, di = 3, n = 60 




Figure 22. Inverse Kasteleyn matrix entries as a function of u> , Wi for 
do = 1, di = 4, n = 60 




Figure 23. Inverse Kasteleyn matrix entries as a function of wq, w% for 
d = 1, di = 4, n = 100 




Figure 24. Inverse Kasteleyn matrix entries as a function of wq, w± for 
do = 2, di = 4, n = 100 




Figure 25. Probability of occurence of the shape in figure |27| as a function 
of position, for the Aztec diamond of side 40. 



:(!> 




Figure 26. Probability of occurence of the shape in figure as a function 
of position, for the Aztec diamond of side 40. 





Figure 27. L-shape Figure 28. L-shape 
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